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ABSTRACT 

Aims. We check how the Fermi/LAT data constrain the physics of hot accretion flows that are most likely present in low-luminosity 
AGNs. 

Methods. Using a precise model of emission from hot flows, we studied the flow y-ray emission resulting from proton-proton inter¬ 
actions. We explored the dependence of the y-ray luminosity on the accretion rate, the black hole spin, the magnetic field strength, the 
electron heating efliciency, and the particle distribution. Then, we compared the hadronic y-ray luminosities predicted by the model 
for several nearby Seyfert 1 galaxies with the results of our analysis of 6.4 years of Fermi/LAT observations of these AGNs. 

Results. In agreement with previous studies, we find a significant y-ray detection in NGC 6814. We were only able to derive upper 
limits for the remaining objects, although we report marginally significant (~ 3 (t) signals at the positions of NGC 4151 and NGC 
4258. The derived upper limits for the flux above 1 GeV allow us to constrain the proton acceleration efliciency in flows with heating 
of electrons dominated by Coulomb interactions, which case is favored by the X-ray spectral properties. In these flows, at most 
~ 10% of the accretion power can be used for a relativistic acceleration of protons. Upper limits for the flux below 1 GeV can 
constrain the magnetic field strength and black hole spin value; we find these constraints for NGC 7213 and NGC 4151. We also 
note that the spectral component above ~ 4 GeV previously found in the Fermi/LAT data of Centaurus A may be due to hadronic 
emission from a flow within the above constraint. We rule out this origin of the y-ray emission for NGC 6814. For models with a 
strong magnetohydrodynamic heating of electrons, the hadronic y-ray fluxes are below the Fermi/LAT sensitivity even for the closest 
AGNs. In these models, nonthermal Compton radiation may dominate in the y-ray range if electrons are efficiently accelerated and 
the acceleration index is hard; for the index == 2, the LAT upper limits constrain the fraction of accretion power used for such an 
acceleration to at most ~ 5%. Finally, we note that the three Seyfert 2 galaxies with high starburst activity NGC 4595, NGC 1068, 
and Circinus show an interesting correlation of their y-ray luminosities with properties of their active nuclei, and we discuss this in 
the context of the hot flow model. 
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1. Introduction 

Low-luminosity AGNs are probably powered by optically 
thin, hot accretion flows (for recent review s see, e.g., 
Poutanen & Veledin^l2014tlYuan & Naravanll2014ll . In the cen¬ 
tral parts of these flows, protons have energies above the thresh¬ 
old for pion production, and then the flows produce consider¬ 
able y-ray fluxes that are due to the decay of neutral pions. This 
property was noted ear ly in the developme nt of the accretion 
theory, for instance, by IShaniro et al.l (Il976h . and was reported 
for the currently ap plied advection-dominat ed accretion flow 
(ADAF) models by iMahadevan et li] JT^. As estimated in 
IShaptro et al.l (Il97^ and lOka & Manmotol (l2003l) . the y-ray lu¬ 
minosity, Ly, of a flow surrounding a rapidly rotating black hole 
may be similar to its X-ray luminosity, Lx- However, these stud¬ 
ies neglected the general relativistic (GR) transfer of y-ray pho¬ 
tons as well as the absorption of these photo ns in the radiation 
field o f the flow. Improved computations by iNiedzwiecki et alJ 
(l20I3l) . hereafter NI3, showed that Ly is at least an order of mag¬ 
nitude lower than Lx if these effects are taken into account. The 
precise y-ray emission model, however, has not been directly 
compared with current observational data. 


No prominent y-ray signal was found based on obser- 
vations in the Fenni/LAJ' data from radio-quiet AGNs (e.g., 
lAckermann et al.l l2012aL l2015ll . with a few notable excep¬ 
tions. These include the three X-ray bright est Seyfert 2 galax¬ 
ies, NGC 4945, NGC 1068, and Circinus (iLenain eF^ 
lAckermann et al.l 12012bl : iHavashida et al.l l2013h . However, all 
three are composite AGN plus starburst systems, and it has not 
yet been definitively established whether their y-ray emission 
is related with the starburst or with the AGN activity. The up- 
pe r limits (UL) on t he pho ton flux in the GeV range, derived 
in lAckermann et al.l (1201 2al) . typically constrain the luminosity 
ratio to L^/Lx < 0.1 and in several Seyferts to Ly/Lx < 0.01, 
which shows that the sensitivity of the Fermi/LAT surveys has 
reached the level at which predictions of the hadronic emission 
from hot accretion flows can be probed. This motivated the de¬ 
tailed comparison we present here. 

We analyzed 6.4 years of Fermi-LAT data of nearby, 
low-luminosity AGNs, which means that we use a data set 
that is more tha n three years longer than the one used by 
lAckermann et al.l (1201 2ah and is more than two years longer 
than the recent third ca talog of AGNs detected by Fermi/LAT 
(lAckermann et al.ll2015l) . A detailed comparison with the model 
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requires a precisely determined black hole mass and intrinsic 
X-ray luminosity, therefore we focus here on several of the 
best-studied objects. As a particularly interesting case, we thor¬ 
oughly examine the data fro m NGC 4151, for which th e ratio 
of LyjLx < 0.0025, found in lAckermann et akl (1201 2al) . is the 
lowest of the 120 Seyfert galaxies considered in their work. 

The physics of particle acceleration and heating in hot flows 
is not well understood. If non-thermal processes (e.g., magnetic 
reconnection) take place, a power-law-like distribution of par¬ 
ticles may be expected. We here take into account predictions 
of hot flow models with proton distributions that include ther¬ 
mal and nonthermal components. We compute the y-ray emis¬ 
sion from accretion flows following the approach of N13. We 
use a GR hydrodynamic description of the flow combined with 
GR, Monte Carlo (MC) computations of radiative processes. We 
improve the model of N13 by refining the hydrodynamic solu¬ 
tions for flows whose proton distribution is dominated by a non¬ 
thermal component. We obtain solutions for a range of accretion 
rates that allows comparing our results with objects that have 
bolometric luminosities of between ~ 10 and 10“^ of the Ed¬ 
dington limit. In Sect. |2 we discuss the dependence of the pre¬ 
dicted Ly on various parameters of the hot flow model, this has 
not been fully explored in precise computations before. In par¬ 
ticular, N13 only considered a weak magnetic held and a fixed 
accretion rate. In Sect. [2 we present the results of our analysis 
of the FermifLAT data together with some results collected from 
the literature. Our main results are presented in Figs.[T]and[2 In 
Sect. |4]we briefly consider effects that are related with the non¬ 
thermal acceleration of electrons, in particular, their Compton 
emission in y-rays. In Sect. |5]we discuss the general implica¬ 
tions of our results for the physics of hot flows; we also briefly 
discuss hard spectral states of black hole binaries, where similar 
hot flows are most likely present. 


2. Model of X- and y-ray emission from hot flows 


We applied the m odel developed in Niedzwiecki et al.l (12012i 
120131120l412015l) and IXie et akl (l201ol) . Here we briefly sum¬ 
marize its basic properties and some improvements we added 
here. 

We considered a black hole that is characterized by its mass, 
M, and angular momentum, J, surrounded by a geometrically 
thick accretion flow with an accretion rate, M. We defined the 
following dimensionless parameters; r - R/Rg, a = Jf(cRgM), 
m - MjM-Eid, where M^dd - R^ddlc^, Rg - GMjc^ is the gravi¬ 
tational radius and LEdd = TnGMm^clerr is the Eddington lumi¬ 
nosity. We assumed a viscosity parameter, a - 0.3 and the ratio 
of the gas pressure (electron and ion) to the magnetic pressure, 

P- 


The fraction of the dissipated energy that directly heats elec¬ 
trons is denoted by 6. We considered two limiting cases, one 
with d = 10'^, for which electron heating is dominated by 
Coulomb interactions at higher m and compression at lower ih, 
and the other with 6 = 0.5, for which electron heating is dom¬ 
inated by magnetohydrodynamic (MHD) processes. Our results 
for 6 - 10“^ are valid for all values of d $ 0.1. 

The Eddington ratio for a given energy range is denoted by A 
with a relevant subscript, for example, A 2-10 kev = L 2-10 kev/Tedd 
and /li_ioGeV = Ti -10 Oev/CEdd, where L 2 -iokeV and Li-iooev 
are the luminosities in the 2-10 keV and 1-10 GeV range de¬ 
tected far away from the flow, that is, they are affected by GR 
transfer and yy absorption effects (only the dashed line in Eig. 
[TJ) shows the rest-frame luminosity). Similarly, the photon flux 



Fig. 1. (a) X-ray photon spectral index as a function of intrinsic 
di-iokev- The observational data correspond to the best fits of the high- 
quality observations described in Sect.|2 The gray regions indicate the 
location of solutions for the model with seed photons from thermal syn¬ 
chrotron only (lower) and including nonthermal synchrotron (of pion- 
decay or directly accelerated electrons; upper panel). Ticks in the top 
axis show the approximate value of Gr (see text). The thick vertical 
line indicates the maximum luminosity for the standard ADAF solu¬ 
tions (i.e., unaffected by the Coulomb cooling of protons), (b) /1 i_iogcV 
as a function of /l 2 _iokcV 5 see Tables[T]and|2for the observational data. 
The solid lines are for model N with 6 = 10“^, a = 0.95, p = 9 (tri¬ 
angles), model N with 6 = 0.5, a = 0, /( = 9 is plotted with circles), 
model Ho.i with 5 = 10“^, a = 0.95, p = 9 with squares, and all three 
are the same as in Fig.[2j3 below. The dashed line with open squares de¬ 
notes the rest-frame 1-10 GeV luminosity in the last model (i.e., Ho.i). 
In panel (b), A 2-10 keV corresponds to the average X-ray luminosity dur¬ 
ing the Fermi/LAT observations (see text). This differs in some objects 
by up to a factor of 2 relative to panel (a). The point for Cyg X-1 shows 
the hard state parameters according to lMalvshev et al.l ll2013h . 


integrated over a given energy range is denoted by F with a rel¬ 
evant subscript. 


We found the global hyd rodynamical solu tions of the GR 
structure equations followingjManmo^ 200011 with minor im¬ 
provements as described in iNiedzwiecki et ~l2012h . From this, 
we found the self-consistent electron temperature distribution, 
7’e(r), for which the electron energy balance when global Comp¬ 
ton cooling is achieved. This involved several iterations between 
the solutions of the electron energy equation and the MC Comp- 
tonization simulations. The procedure involves the assumption 
that the flow structure is unaffected by changes in Eg, which con¬ 
strains the range of exact solutions to /l 2 -iokeV $ (1 - 2) x 10“^ 
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(depending on the model, see also I Yang et al.l (l2015h '): the max¬ 
imum /l 2 -iokeV is denoted by /lADAF.max and is indicated by the 
thick vertical line in Fig. [T^. The Coulomb cooling of protons 
becomes strong at higher luminosities, and the flow is char¬ 
acterized by a dramatic dependence on even small changes of 
Te, (when globa l Compton cooling is taken into account, cf. 
IXie et al] (l201()l) '). which should lead to a collapse of the flow 
and formation of a standard cold disk at r > r^, with the transi¬ 
tion radius, rtr, dependent on m. 

All model results we present here were obtained for flows 
extending out to r = 10'*. We currently investigate the regime 
of A 2 - 1 Q kev > dADAF.max- Although the results are rather prelimi¬ 
nary, we And that (1) up to /I 2-10 keV - 3 x 10“^ (indicated by the 
boundary of gray area in Fig.[T^) the hot flow can occupy the in¬ 
ner region within rtr - 100 (this accretion flow geometry is sup¬ 
ported by observations, as noted below); (2) up to this /l 2 -iokeV, 
the X-ray and the rest-frame y-ray emission is correctly approx¬ 
imated by models with flows extending to large r, for which we 
used the hydrodynamical solution neglecting the Coulomb cool¬ 
ing of protons; and (3) this model (with a flow extending to large 
r) may overestimate the y-ray opacity (see below). 

The minimum luminosity of flows that can be studied with 
our computational method, determining the lowest /I 2-10 keV (de¬ 
noted by dniin) in all presented models, is set by the requirement 
that the energy balance for electrons is determined by radiative 
cooling (rather than advection). Close to Tnjin , the compressive 
heating of electrons dominates for 6 = 10“^, yielding an approx¬ 
imately constant radiative efficiency. This effect is crucial for as¬ 
sessing the y-ray emission in objects with A 2 - 1 Q kev ^ , and 

we return to this below. 

The magnitude of the hadronic processes involves a signifi¬ 
cant uncertainty related to the distribution of proton energies. At 
d 2 -iokeV < dADAF.max the thermalization and cooling timescales 
for protons are much longer than the accretion timescale, and 
the distribution of protons is determined by poorly understood 
heating and acceleration processes. We considered three cases: 
a thermal model T in which all protons have a Maxwellian dis¬ 
tribution; 

a nonthermal model N, where the total power is used to acceler¬ 
ate a small fraction of protons, for which we assume a power-law 
distribution, «pi(y) oc y~^ up to y^ax = 100, and we assume that 
the other protons remain cold; the radius-dependent fraction of 
the protons with the power-law distribution is given by Eq. (7) 
inN13; 

and a hybrid model Hq.i that has 10% of the energy content in 
nonthermal protons, denoted by T/p, and 90% of the energy con¬ 
tent in thermal protons; for the nonthermal protons we assume a 
power-law distribution like that in model N. 

The assumptions of these models are identical to models T, N, 
and H in N13, except for the energy content of nonthermal com¬ 
ponent in model H, which was radius-dependent in N13, whereas 
here we assumed a radius-independent rjj, = 0.1 (indicated by 
the subscript in Ho.i). For model N, however, we here take into 
account the effect neglected in N13, namely, a higher compress¬ 
ibility of the flow supported by the pressure of a nonthermal, 
relativistic proton gas, corresponding to the decrease of the adi¬ 
abatic index from ^5/3 (for a thermal proton gas with typically 
subrelativistic temperatures) to 4/3. All results for model N in 
this work use the hydrodynamic solutions including this effect. 
The higher compressibility of the flow results in its smaller geo- 
metrical thickness and, hen ce, higher density (see also Fig. 4 in 
iPopham & Gammii^ (Il998h ). For a given set of parameters, the 
density is higher by a factor of ~ 1.5 (for /3 - 1) to ~ 5 (for 
jS = 9) in model N than in models T or Ho.i. 


We note that the dependence of the ge ometrical structure 
on P in model N, which was discussed in iNiedzwiecki et al.l 
Jml . is almost completely reduced. The results presented for 
this model, with both components significantly contributing to 
the total pressure (i.e., protons and magnetic held), described by 
the adiabatic index 4/3, correspond to the maximum density, and 
hence the maximum Thomson depth, t, as well as the maximum 
Ly of the flow for a given set of ih, a, jA, and 5. 

To compute the X-ray radiation, we took into account seed 
photons from nonthermal synchrotron emission, which is crucial 
to reach agreement for the model spectra with th e observational 
data. As first noted by lYuan & Zdziarski (120041), ADA F mod¬ 
els that follow the original formulation by Naravan & Yil (Il99-5h . 
that is, whose seed photons are provided by thermal synchrotron 
emission, predict too high and disagree with observational 
data f or AGNs. Following the results of INiedzwiecki et al.l 
dMl, we self-consistently included the synchrotron emission 
of pion-decay electrons for 6 - 10“^, whereas for S - 0.5 we as¬ 
sumed that a fraction of the electron heating power was used 
for their nonthermal acceleration. In the latter case, we assumed 
//e = 0.1 and a monoenergetic injection of electrons at y = 100, 
for which synchrotron self-absorption is relatively unimportant 
(in Sect. 14.21 we release this assumption and consider the depen¬ 
dence on the electron acceleration distribution). The two gray re¬ 
gions are delineated in Fig. [T^ to indicate the location of model 
solutions that include a significant nonthermal synchrotron radi¬ 
ation (upper panel) and solutions with only thermal synchrotron 
emission (lower panel), see Sect. |5] for a more detailed discus¬ 
sion. 

Using our MC Comptonization code, we tabulated the 
angular-, energy-, and location-dependent distribution of pho¬ 
tons propagating in the central region. Then, to compute Ly , 
we took into account the absorption of escaping y-ray photons 
on pair creation in interactions with the (tabulated) target pho¬ 
tons, strictly following N13. In Fig. |2^ we show the size, rpj,, 
of the y-ray photosphere (inside which the flow is opaque to y- 
rays), which is determined as the radius of the emission point in 
the equatorial plane, from which the optical depth along the out¬ 
ward radial trajectory is Tyy - 1. For radial directions close to the 
symmetry axis, values of rph determined by a similar condition 
are lower by a factor of ~ 1.5. Results presented in this work, 
except for the dashed line in Fig. [2 correspond to spectra aver¬ 
aged over the observation angle, flobs- A strong dependence on 
0obs only occurs for a rapidly rotating black hole and only at low 
42-iokeV, at which the contribution from the innermost few Rg 
is not attenuated by yy absorption. In other cases, the difference 
between y-ray fluxes received at different 0obs does not exceed a 
factor of ~ 2. The amount of attenuation due to yy absorption is 
illustrated by the dashed (rest-frame Ly) and solid (observed Ly) 
lines with squares in Fig.[TJ); we note that at high /l 2 -iokeV the 
internal Ly is higher by two orders of magnitude than observed. 

The dependence of rph on /l 2 -iokeV is similar in all models 
(i.e., for different a, /3, or 6). However, the scaling presented 
in Fig. |2^ cannot be simply extrapolated to /l 2 -iokeV ^ 10“^, 
at which this extrapolated rph would exceed the flow size of 
rp - 100. In this case, the y-ray opacity should be compara¬ 
tively lower, at least in outer parts of the flow (i.e., close to rp). 
In a flow extending to large r, the opacity for a y-ray photon 
escaping from rem is dominated by head-head interactions with 
X-ray photons produced at r > repj. For photons escaping from 
the outer parts of a flow contained within rp, the opacity is dom¬ 
inated by much less efficient head-tail interactions with X-ray 
photons produced at r < rem. 
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Fig. 2. (a) Photosphere radius for photons with £ = 0.1, 1, and 10 GeV from bottom to top as a function of /l 2 _iokeV in the model with 6 = 10“^, 
a = 0.95, y6 = 1 (in other models the dependence of rph on /l 2 -iokeV is very similar); (b) /li_ioGeV as a function of /i 2 -iokcV in models with r/p > 0 
and i = 2.1; the dotted lines plot models N with d = 10“^, (magenta) circles indicate models N with a = 0.95 and/3 = 9 with m = 0.05, 0.1, and 
0.3, (green) squares stand for a = 0, = 9, m = 0.1, 0.3 and 0.5; the dashed lines show models Ho.i with 6 = 10~^, (black) circles models with 

a = 0.95, = 1 with m = 0.033, 0.1, 0.3 and 0.5; (blue) squares indicate a = 0.95, /? = 1 with m = 0.1, 0.3, 0.5, 0.8 and 1.2; the solid lines show 

models N with d = 0.5, (red) squares depict a = 0, /? = 1 with m = 0.004, 0.01, 0.033, 0.1 and 0.3; (orange) circles stand for a = 0.95, /? = 9 with 
m = 0.003, 0.01, 0.033 and 0.1. (c) The high-energy part of the observed spectra of leptonic emission. (Black) solid lines indicate model Ho.i with 
S = 10“^, a = 0.95,/3 = 1 and m = 0.033, 0.1, 0.3; (red) dashed lines show model N with 6 = 0.5, a = 0,/3 = 1 and m = 0.004, 0.01, 0.033, 0.1; 
(blue) dotted line shows model N with 6 = 0.5, a = 0.95,= 9 and m = 0.01. 



Fig. 3. Ao.2-1 Gev as a function of /t 2 _iokev> see Tables [T] and |2] for the 
observational data. All model points are for thermal protons (model T) 
with 6 = 10“^; the solid line with circles shows a = 0.95, /3 = 1, 
lit = 0.033, 0.1, 0.3, and 0.5; the solid line with full squares depicts 
a = 0.95, p = 9, m = 0.1, 0.3, 0.5, 0.8, and 1.2; the dashed line with 
open squares plots the estimate (see text) for an edge-on observer for 
a = 0.998, yS = 9, m = 0.033, and 0.1. 


Below we summarize the dependence of Ly on various pa¬ 
rameters that is relevant for the comparison presented in Sect. [3] 
For thermal protons, the bulk of the hadronic radiation is emit¬ 
ted in the 0.1-1 GeV range. The presence of nonthermal protons 
is reflected in the production of photons with £ > 1 GeV, with 
a power-law spectrum for a power-law proton distribution. We 
consider these two regimes separately. 

2.1. Thermal protons: Ao. 2-1 cev vs 6, p, and a 

The rate of pion production depends on the number of protons 
with energies above the pion production threshold. For the ther¬ 
mal distribution of protons, this rate is extremely sensitive to 
the proton temperature, Tp. The pion production rate is low for 
Tp < 10*2 K and negligible for Tp < 4 x 10*' K. The highest 
Tp, and accordingly the highest Ly, can be expected for flows 
with large p and a and small 6. The maximum /I 0 . 2-1 GeV that 


can be observed for thermal protons, with /I 0 . 2-1 GeV - 10“^ for 
T 2 -iokeV between ~ lO^'* and 10 “^^ is shown by the line with 
squares in Fig. [3 In models with a = 0, 6 - 0.5 or p = 1, 
do.2-1 GeV is smaller at least by a factor of several, and even by 
orders of magnitude at /i 2 -iokeV - 10 “-*; in all these cases, the 
model predictions are below the Fermi ULs discussed in Sect. [3] 

The dependence on these three parameters involves the fol¬ 
lowing effects. In flows with smaller p, a larger part of the ac¬ 
cretion power is used to build up the magnetic held strength, 
therefore, the proton heating power is weaker. As a result, for 
P - I the proton temperature is by a factor of ~ 2 lower than for 
P - 9. The related decrease of /I 0 . 2-1 GeV is illustrated in Fig. [3 
for = 1 and a - 0.95, Tp only exceeds 4 x 10** K at r < 6 , 
hence, at /l 2 -iokeV > 10 “"* the y-ray emission is strongly atten¬ 
uated by 77 absorption. The increase of b yields the increase of 
electron heating at the expense of the proton heating, and it there¬ 
fore obviously results in the decrease of the £y/£x ratio. For the 
decreasing a, a decreasing fraction of the accreted rest-energy is 
dissipated, therefore, Tp decreases. In the extreme case of a = 0 
and p - \,Tp weakly exceeds 4 x 10** K only close to the event 
horizon, and no hadronic emission can be expected. 

At 42-10 kev - 10“^ photons with E % \ GeV may escape 
from the ergosphere (the extremely relativistic region located 
within r < 2) and some specific effects of the Kerr metric can be 
observed. These are that gravitational focusing of the photon tra¬ 
jectories toward the equatorial plane enhances the flux received 
by edge-on observers. Furthermore, the increase of a from 0.95 
to 0.998 increases Ly by a factor of ~ 2; at r > 2 both values 
of a give the same 7 -ray production rate. Taking these effects 
into account, we show the estimated maximum Ao. 2-1 GeV that 
can be expected at /l 2 -iokeV ~ 10“^ by the dashed line in Fig. 
13 For p - 9 this /l 2 -iokeV is lower than Tmin, meaning that it 
corresponds to the regime where advection of electron energies 
dominates. This in turn means that the lowest /I 2-10 kev point (for 
a - 0.998 and m - 0.033) is not based on a precise solution of 
our mod el; we used the sca ling of radiative efficiency, oc m^ '’, 
found in IXie & Yu^ (1201 2h for low m, to roughly estimate the 
X-ray luminosity. 
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2.2. Nonthermal protons: /li-iocev vs 6 and t/p 

The dependence of /li_io cev on A 2-10 kev in models with t/p > 0 
is shown in Fig. [JJt. Here the number of protons above the pion 
production threshold depends only linearly on the proton heating 
power. This means that Ly is much less sensitive to both (3 and a 
than in model T. The difference between various models results 
mostly from different radiative efficiencies (of electrons), which 
establishes the relation between /I 2-10 keV and m, and accordingly 
the rest-frame Ly oc n-p-. The efficiency is determined primarily 
by the value of 6. Obviously, /ii_io ceV also depends on //p. How¬ 
ever, we note that the difference between models Hq.i and N does 
not simply involve the difference of /li_ioGeV by a factor of 10, 
corresponding to an increase of from 0.1 to 1. The difference 
is more complex because of the difference between the hydrody- 
namical solutions that underlie these models. 

For 6 - 10“^, when the heating of electrons is dominated 
by Coulomb interactions (at larger m), the rest-frame y-ray lu¬ 
minosity scales roughly linearly with the X-ray luminosity (both 
oc m^). For this scaling, combined with the increase of rph with 
increasing m, the observed /li_ioGeV is roughly constant while 
d 2 -iokeV changes by more than an order of magnitude (see dot¬ 
ted and dashes lines in Fig.lJJt). At low m, when the compressive 
heating of electrons dominates, leading to /I 2-10 keV we note 
a fast decrease of /li_io ceV with decreasing A2-10 kev; this is most 
clearly seen for model Ho.i with (3-9 (dashed line with squares 
in Fig.|2j3), in which compressive heating dominates for m % 0.3, 
corresponding to/i 2-10 keV ^ 10 ^“*. 

In models with 6 - 0.5, the radiative efficiencies are higher, 
therefore, /ii_io Gev is lower (because of the lower m) at a given 
d 2 -iokeV than for 6 - 10^^. Hence, for 6 = 0.5 we only present 
results for model N; for Hq 1 the predicted y-ray flux is below 
the sensitivity of LAT. For a constant 6 , A2-10 kev results in 
a fast decrease of /li_ioGeV with decreasing /l 2 -iokeV- Further¬ 
more, for large 6, the radiative efficiency increases with a, as a 
result, GeV at a given /I 2-10 keV is by a factor of several larger 
for a - 0 than for a - 0.95 (somewhat counter-intuitively). 

The uncertain value of 6 introduces the main uncertainty in 
testing the hadronic emission model by allowing for a range of 
m at a given /l 2 -iokeV- However, different m also correspond to 
different Thomson depths of the flow, which is reflected in the 
shape of the X-ray spectrum. The above uncertainty can there¬ 
fore be reduced by using the X-ray spectral information. Figure 
12}; shows example X-ray spectra for large and small 6. Spec¬ 
tra of large-(5 models in general deviate from a power-law shape 
at /i 2 -iokeV ^ 10“"^, that is, for m ^ 0.01 with a - 0.95 and 
m ^ 0.03 with a = 0. For low-(5 models, deviations from a 
power-law shape are weak down to /l 2 -iokeV ~ 10“^ (the bot¬ 
tom solid line in Fig.|2};). 

We considered various proton distribution slopes between 
s = 2.1 and 2.7. In all models, the lowest di-ioGeV corresponds 
to i = 2.1 (for s - 2.6, /li_ioGeV is larger by up to a factor 2), 
and conservatively we assumed this s for Fig. [It. 


3. Sample and data analysis 

We considered nearby AGNs with spectral properties consistent 
with the hot flow scenario. This includes several Seyfert 1 galax¬ 
ies: NGC 4151, NGC 5548, IC 4329a, NGC 6814, NGC 4258, 
and NGC 7213. We also took into account the FR I galaxy Cen- 
taurus A, whose X-ray radiation may be dominated by a Seyfert- 
like emission. Spectral parameters of the intrinsic X-ray emis¬ 
sion, measured in all these objects, agree with our model of ther¬ 
mal Comptonization in hot flows (see, e.g., FigUt). Furthermore, 


Table 1. Adopted 2-10 keV Eddington ratio, black hole mass, distance, 
and references for the data. The references are given in the order X-ray 
measurement, M, and (if available) distance. 


Source 

42-10 keV 

xlO^^ 

M 

lO^Mo 

D 

(Mpc) 

Refs. 

IC 4329a 

0.48 

13™ 

68.4 

ffi,2 

NGC 4151 

0.12 

5.4+ 1.8 

19.0 

3* ,4,4 

NGC 4258 

0.00094 

3.6 

7.2 

5" ,6,7 

NGC 5548 
NGC 6814 

0.28 

0.15 

3 9+2-3 

026-8-.^io 

72.7 

22.8 

8,9 

10"',9,11 

NGC 6814* 

0.021 

1.9 ± 0.4 

22.8 

10"',12,11 

NGC 7213 

0.014 

8 ■0-^0° 

22.0 

13,14,11 

Cent. A 

0.012 

5.5 ±3.0 

3.8 

15,16,17 

Circinus 

1.1-2.3 

0.17 + 0.03 

4.2 

18,19,11 

NGC 1068 

1.7 

1.0 

14.4 

20,21,11 

NGC 4945 

1.8 

0.14 

3.6 

22,23,24 


Notes. A 2 -W kev given in the second column involves rescaling (based on 
the 5w!/f/BAT data, see text) of the average intrinsic luminosity given in 
the reference by a factor of 1.3; ® 1.2; 2; 1.5. 


References. (1) 

Brenneman et all (I2014h. 121 iMarkowitzl (120091). 

131 iLubihski et al. 

(201011, 141 

Hbnie et all (120141). 151 Revnolds et al. 

(I 2 OO 9 I). 161 Mivoshi et al.l 

19951). 171 iHermstein et all i 

1999b, 

(8) 

iBrenneman et all (120121). 191 iPancoast et al. i 

201.5h. 

101 

Walton et all (l2013ll. 1111 iTullvl (Il988ll. 1121 iBentz et al. 

200? 

). 1131 Lobban et al. 

(|201Qh. 1141 ISchnorr-Muller et al. 


inUHarris et allll2010h. ('18 ) lArevalo et'^feOldli. tlQllGreenhi 1 et all 

( I2003t). (201 ' ~lBauer et al. I ( I2014li. (21) iGreenhill et all (Il996t). 
' Puccetti et all(l2014li . r23i lGreenhill et alldlQQTh . t24i lTullv et all 


they show direct evidence of the lack of an optically thick disk 
in the central region. For objects with the lowest luminosities 
(NGC 4258 and NGC 7213) studies of the Fe Ka line place 
the in ner edge of optically thick disk between rtr ~ 10^ and 
~ 10^ (iRevnolds et al.ll2()0^ iLobban et al.ll201(lh . In more lumi¬ 
nous objects (NGC 4151, NGC 5548, IC 43 29a, and NGC 6814), 


the measured Ee Ko; line widths, % 100 eV ( 

Lubinski el al.l2010l: 

IWalton et al.ll2013l IBrenneman et al.ll2012l 

2014b. indicate trun- 


cation of the disk at rn ^ 100; a similar location of an op- 
tically thick disk is indicated by reverberation measurements 
(lEdelson et al.ll20l'^ . These accretion geometries are consistent 
with the model of spect ral evolution deve loped for black hole 
binaries, for instance, bv lEsin et al.l (1 19971). with r,r changin g in 
response to changes in m (for a review see lDone et al.ll20(y^ . 

We also considered three Seyfert 2 galaxies, NGC 1068, 
Circinus, and NGC 4945. All three are observed to have high 
bolometric luminosities, above the range that can be precisely 
studied with the current version of our model. Nevertheless, 
comparing their y-ray luminosities with properties of their ac¬ 
tive nuclei appears to be interesting. 

Eor most AGNs in our sample, direct measurements of the 
masses of their supermassive black holes are available, only for 
IC 4329a and NGC 7213 did we use estimates based on the ve¬ 
locity dispersion. The adopted M and distance values are given 
in Table |T] Eor NGC 5548 and IC 4329a we assumed the lu¬ 
minosity distances for the following cosmological parameters: 
Gm = 0.27, Ga = 0.73, and Hq - 1\ km s ' Mpc *. The error 
bars on the values of A in Eigs. [T]and|3]are mostly due to uncer¬ 
tainties in the determination of M. Of the measurements of M 
that are based on H 2 O megamaser kinematics, the uncertainty is 
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Table 2. Integrated photon flux, F, and the photon index, F, of power- 
law fits, or 95% confidence level ULs with the assumed F, and the cor¬ 
responding Eddington ratio for 1-10 GeV (upper part) and 0.2-1 GeV 
(lower part). 


Source 

F 

F 

/iy 


10 

“'(* ph/cm2/s 

xlO 


Energy range: 1- 

-10 GeV 


IC 4329a 

2.7 

< 0.46 

< 0.49 

IC 4329a 

2.1 

< 0.68 

< 0.88 

NGC 4151“ 

2.7 

< 1.2 

< 0.23 

NGC 4151“ 

2.1 

< 1.3 

< 0.31 

NGC 4151* 

2.1 

< 1.8 

< 0.43 

NGC 4151“ 

2.1 

< 1.8 

< 0.43 

NGC 4151“' 

2.1 

< 2.0 

< 0.48 

NGC 4258 

2.7 

< 1.1 

< 0.043 

NGC 4258 

2.1 

< 1.5 

< 0.067 

NGC 5548 

2.1 

< 0.3 

< 1.8 

NGC 5548 

2.7 

< 0.6 

< 2.9 

NGC 6814 

2.6+ 0.1 

3.0+ 0.7 

18 + 10 

NGC 6814* 

2.6+ 0.1 

3.0+ 0.7 

2.6+1 

NGC 7213 

2.7 

< 0.4 

< 0.07 

NGC 7213 

2.1 

< 0.5 

< 0.11 

Cent. A^^' 

2.1 +0.2 

18.0 + 2.7 

0.17!°-2 

Circinus^^^ 

2.4+ 0.1 

9.2 + 1.4 

3.1 +0.7 

NGC 1068(2) 

2.3+ 0.1 

5.4+ 0.8 

3.7 + 0.8 

NGC 4945(2) 

2.4+ 0.1 

9.8 + 1.2 

2.9 + 0.6 

Energy range: 0.2-1 GeV 

IC 4329a 

4.0 

< 3.3 

< 0.51 

NGC 4151 (1)“ 

4.0 

< 4.5 

< 0.13 

NGC 4151 (2)* 

4.0 

< 20 

< 0.52 

NGC 4258 

4.0 

< 7.7 

< 0.044 

NGC 5548 

4.0 

< 12 

< 8.4 

NGC 7213 

4.0 

< 7.1 

< 0.18 


Notes. using the full data set of 6.4 years and including source S in 
the model; ~ 4.9-year data set (6.4 years without period S), source S 
not included; data from the first three years, S included in the model; 
data from the first three years, S not included in the model. 

Re ferences. _ 

(D lSahak yan et al.](|2013h. parameters of the ’second’ component above 
4 GeV: (2l lAcero et alJll20lli . 


reported only for Circinus; for the remaining cases (NGC 1068, 
NGC 4258, and NGC 4945) we assumed the same relative un¬ 
certainties in M for Figs. [T]and[3] Fo r NGC 6814 the d iffer¬ 
ence of measuremen ts of M between iBentz et ^ (l2009h and 
iPancoast et al.l (1201 5h amount to an order of magnitude, so we 
present results for both values; the case with the higher value of 
M is denoted by a star superscript in Tables [1] and |2 and Fig.[T] 

Table |2] shows the results of our analysis of the FermifLAfT 
data, described below, except for NGC 1068, Circinus, NGC 
4945, and Centaurus A, for which we used the results of previ¬ 
ous studies collected from literature. In Figs.[TJ) and|3]we show 
the measured /I 0 . 2-1 cev and /li_io cev, or UL on them, as a func¬ 
tion of the intrinsic A 2-10 kev- To find the latter, we used results of 
detailed X-ray spectral studies (specified below) of high-quality 
data from Suzaku, NuStar, or simultaneous XMM and INTE¬ 
GRAL observations, which allow distinguishing the primary X- 
ray emission and reflection or absorption components. The pa¬ 
rameters of these fits are shown in Fig.[T^. Using the 5w!/f/BAT 



Time [MJD] 


Fig. 4. Small black circles with error bars show SwiftfBKf 15-50 keV 
light curves in seven-day bins for NGC 7213, NGC 4151, IC 4329a, and 
NGC 6814. The red vertical line in each panel shows the average flux 
during the 6.4 years of analyzed FermifLKf observations. The large 
orange circles in panels (acd) show the times of the adopted Suzaku 
and NuStar observations and the average flux in time bins equal to the 
duration of these pointing observations. The red vertical line in panel 
(b) indicates the end of th e three-year period used for comparison with 
lAckermann et all (l2012ah . The blue line in the bottom of panel (b) de¬ 
lineates period S (see Sect. 13.21 . Th e flux levels of dim ( D), medium 
(M), and bright (B) states, defined in lLubinski et all ll2010h . are shown 
in the left part of panel (b). 

(iKrimm et al.l2013h light curved , we find that even for the most 
variable objects (NGC 7213, IC 4329a, and NGC 6814; their 
light curves are shown in Fig. H the adopted X-ray observa¬ 
tions correspond to the flux levels close to the average 5wt/f/BAT 
flux during the analyzed 6.4 years of FermifLPff observations; in 
some objects we used the intrinsic 42-10 keV rescaled to account 
for a small difference, see Table [T] 

3.1. LAT data analysis 

For NGC 4151, NGC 5548, NGC 4258, NGC 6814, NGC 
7213, and IC 4329a we analyzed 6.4 years of the Femti/LAT 
data, comprising observations carried out between 2008 Au¬ 
gust 4 and 2015 January 10. For each object, events were se¬ 
lected from a region with a radius of 15° centered on the posi¬ 
tion of the analyzed source. We performed the unbinned like¬ 
lihood analysis using the v9r33p0 Fermi Science Tools with 
CALDB instrument response functions. We used the standard 
templates for the Galactic (gll_iem_vQ5_revl. fits) and the 
isotropic (iso_source_vG)5_revl. txt) backgrounds. In the 
likelihood analysis we took into accou nt all sources repor ted in 
the FermifLAfT Third Source Catalog (lAcero et alJl2015i here- 


* http://swift.gsfc.nasa.gov/results/transients 
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after 3FGL) within a radius of 15° around the analyzed object. 
Each of the catalog sources was modeled with a best-fit spectral 
function (as specified in the catalog, in most cases a power-law) 
with parameters left free in the model fitting. 

Except for NGC 6814 we did not find statistically significant 
signals. We derived the 95% confidence level ULs for the photon 
flux using the Bayesian method. 

We first checked how the use of both an extended data set and 
improved models (based on 3EGL) of the celestial regions of in¬ 
terest affects the UL values. Assuming the y-ray photon spectral 
index E = 2.5, we found ULs for E>o.i cev lower by a factor ~ 4 
(for NGC 7213 and I C 4329a) and ~ 2 ( for NG C 5548) than the 
corresponding UL in I Ackermann et all (l2012ah . Eor NGC 4151 
the UL value is somewhat dependent on the approach to model¬ 
ing (as discussed below), but in general, we were unable to re¬ 
duce it significantly below the value quoted by I Ackermann et ^ 
(l2012ah . 

The pion decay spectra can only be described by a simple 
power-law in limited energy ranges. Lor the thermal distribution 
of protons, the jr'^-decay spectrum can be very roughly approxi¬ 
mated by a very soft power-law, with L ^ 4, above ^ 0.2 GeV. 
Lor the power-law proton distribution, it is a power-law above 
~ 1 GeV. Then, to compare with predictions of model T, we as¬ 
sumed L = 4 and found UL for Lo. 2-1 CeV- To compare with pre¬ 
dictions of models with > 0, we found UL for Li_ioGeV as¬ 
suming different values of L between 2.1 and 2.7. Lor all AGNs, 
except for NGC 5548, the highest ULs correspond to L = 2.1 
(the difference for other values of L does not exceed a factor 
of 2), and we used them in Lig. [TJt. The comparison in Lig. [TJt 
is presented in the most conservative manner, meaning that we 
used the lowest /li_ioGeV predicted by the model (for s - 2.1) 
and the highest UL value (for L = 2.1). 

3.2. NGC 4151 

At the position of NGC 4151 we found a signal with the test- 
statistic significance TS = 17 for E - 0.2 - 6 GeV, see Lig.|5}l. 
It most likely contains a contribution from a new y-ray source 
found in our analysis, denoted below by S, shifted by only ~ 0?5 
from NGC 4151. Interestingly, however, the gtfindsrc best- 
fit location for this signal (shown by the black cross), assuming 
that it comes from a single source, is shifted from the position of 
NGC 4151 by only 8'. Source S is strongly variable, which hin¬ 
ders a proper assessment of its contribution to the y-ray signal. 
The details of our analysis are as follows. 

We first examined the significances of y-ray signals around 
NGC 4151 by means of their TS values. This revealed three new 
sources in the region, which were not reported in the FemiifLA^ 
Second Source Catalog a n d henc e not included in the model 
used in lAckermann et al.l (l2012ah . Lor each of these sources 
we used the gtlike and gtfindsrc tools to find their signif¬ 
icance, best-fit position, and spectral parameters. Two of them 
(indicated by the red symbols in Lig.|5^) have recently been re¬ 
ported in 3LGL, with parameters approximately consistent with 
those estimated in our analysis. They are 3LGL J1220.2-1-3434, 
shifted from NGC 4151 by ~ 5°, for which our analysis gives 
TS = 164, r = 2.2 + 0.1 and L>o.iGeV ^ (8.4 + 1.5) x 10^® 
ph s ' cm- 2 , and 3LGL J1203.2-H3847, shifted from NGC 4151 
by ~ 1.5°, for which we derived TS = 32, T = 2.3 + 0.2 and 
L>o.i GeV - (3.5 + 1.4) X 10“® ph s“' cm“^. These sources are 
included in the model used for all results presented in this work. 

The position of the third source, S, ~ 0?5 from NGC 4151, 
is determined by several photons with energies between 10 and 
20 GeV, which arrived from the same direction (within ~8') be¬ 


tween December 2011 and June 2013; we assume that the source 
exhibited an outburst activity during this 1.5-year period, which 
we denote as period S. Source S is clearly seen in the TS map 
built for period S alone, see Lig. whereas the map built with¬ 
out this period only shows a weak residual, centered on NGC 
4151 and not on S, see Lig. |5^. After subtracting all 3LGL 
sources, we find for source S (1) TS = 30, T = 1.78 + 0.37, 
L>o.i GeV - 1-2 X 10“® ph s“' cm“^ using the data from period 
S, and (2) TS = 22, L = 2.17 + 0.20, L>o.i GeV ^ 2.1 x 10^® 
ph s ' cm“^ using the full 6.4-year data set. Source S is not re¬ 
ported in 3LGL, which was built using data up to 2012 July, that 
i, covering only ~ 30% of period S. Lor the data from period 
S, gtfindsrc gives the location of S at cr 72000 = 12''11™27®, 
b 72000 = 38°56'48". A possible candidate for this source is a BL 
Lac object, 2E 1209 .Oh- 3917, only = 4' from this location. 

Lor a 4.9-year data set, without the period S, gtlike shows 
no signal (TS = 1.5) at the S position. The residual seen in the TS 
map for this data set, Lig. |3a), can be fully compensated for by 
adding the point source at the position of NGC 4151 (compare 
with panel b), and gtlike gives TS = 8, L = 2.7 + 0.3 and 
L>o.i GeV = (1.5 + 0.6) X 10“® for this source. 

The presence of a variable source S results in an ambiguity 
for the UL on the y-ray flux from NGC 4151, especially below 
1 GeV, where the point spread function of LAT is much higher 
than the angular separation between S and NGC 4151. In our 
analysis we considered two variants: 

(1) assuming that S strongly softened before and after period S, 
but still provided some contribution at low energies, we used 
the full 6.4-year data set and included source S in the model 
in addition to the 3LGL sources, which yielded the UL values 
denoted by superscript ’a’ in Table |2l 

(2) assuming that S completely faded away before and after pe¬ 
riod S, we used the data for 4.9 years without period S and in¬ 
cluded only the 3LGL sources in the model (i.e., without source 
S), which yielded the UL values denoted by superscript ’b’ in 
Table |2] 

The thus obtained ULs on Lo. 2-1 GeV differ by a factor of ~ 4, 
both values of the resulting UL on Ao. 2-1 GeV are shown in Lig. 
|3] The UL on Li_io GeV is much less model-dependent. In Lig. 
[It we use the UL obtained for the above case (1) with assumed 
L = 2.1, other ULs for this energy range differ by at most 50% 
and do not affect our conclusions. We also note that adding 
source S in the model for the first three years has a negligible 
effect, see cases ’cd’ in Table |2] We also checked that our results 
are unchanged if we introduce additional sources to account for 
several residuals located at ~ 1° to ~ 4° from NGC 4151, seen 
in Lig. Ha). 

We used the recent dust-parallax distance measurement and 
the implied stellar-velocity-based mass (see Table[TJ, both higher 
by ~ 50% than assumed in previous works on high-energy emis¬ 
sion from NGC 4151. NGC 4151 shows a moderate X-ray vari¬ 
ability with changes of the X-ray flux by up to a factor of ~ 4 
between the dim and bright states. Usin g the adopted distanc e 
and M values and the power-law fits of Lubinski et al.l (1201 Oli . 
we found the intrinsic /l 2 -iokeV - (0.5,1,2) x 10“^ for the (dim, 
medium, and bright) state; the parameters of these three states 
are shown in Lig. [T^). The average Swift/BAI flux during the 
considered Leniu'/LAT observations is higher by a factor 1.2 than 
the flux of the medium state, see Lig. St. 

As we see in Ligs [TJ) and |2 the ULs on Lo. 2-1 GeV and 
Li-ioGeV derived for NGC 4151 yield the tightest constraints 
on both Aq, 2 -i GeV and /ii_ioGeV of the bright Seyfert 1 galaxies 
(although IC 4329a has the lowest Ly/Lx ratio in our sample). 
Lor 6 - 10“^, model N overpredicts di-io GeV by a factor of sev- 
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Fig. 5. TS maps for the region around NGC 4151; in all panels the green cross shows the location of NGC 4151 and the orange cross the location of 
source S. In all panels only the 3FGL sources were subtracted from the maps, except for panel (b), where an additional point source was added at 
the position of NGC 4151. (a) The 8° x 8 ° degrees region with a pixel size of 0?2. The map was built for the energy range 0.2 - 30 GeV, neglecting 
period S (see text); the red and white crosses indicate the location of 3FGL sources, (b) The same as in (a), but the model includes a point source 
at the position of NGC 4151 to compensate for the residual seen in (a), (c) The 2?55 x 2?55 degree region with a pixel size of 0?075 for 0.2-30 
GeV for period S. (d) The 4?8 x 4?8 region with a pixel size of 0?075 for 0.2 - 6 GeV and the full period of 6.4 years; the black cross shows the 
location of a point source found with gtfindsrc under the assumption that the residual seen in the map is produced by a single source; source S 
is not seen in this map because photons with £ > 10 GeV are excluded. 


eral, regardless of a or /? values. Then, if heating of electrons is 
dominated by Coulomb interactions, the energy content of non- 
thermal protons is limited to rjp ^ 0.2. 

For 6 - 0.5, the model predictions, except for model N 
with a = 0, are below the UL value. However, as noted be¬ 
fore, the much lower Ly predicted by large-b models is accom¬ 
panied by much lower t. The X-ray spectrum of NGC 4151 
is precisely measured, at least in the bright state, and allows 
testing this property. Fits of the bright-state spectrum with the 
slab corapPS model in iLubihski et all (1201 Ol) give ^ 1.3. To 
compare it with our model predictions, we simulated the X-ray 
spectra for our models with 42-10 kev - 0.002; we used the IN- 
TEGRAL/lSGRl response function and assumed the normaliza¬ 
tion corresponding to the bright state of NGC 4151. Then, the 
simulated spectra were fitted with the slab compPS model. For 
b = 10“^, the best fit has ^ 1.1 (model Hg i with a - 0.95, 


f] — I, m - 0.65). For 6 - 0.5, all simulated spectra were fitted 
with 5 0.5. We conclude that although large-b models are 
essentially unconstrained by the Fermi UL, they also appear to 
be inconsistent with the X-ray data. 


4o.2-i GeV predicted by model T, compared with with the 
Fermi UL for the above case (1), disfavors models with both 
a large a and a large [3, in a manner consistent with NGC 7213. 
For case (2), however, the UL is above the model prediction. 
Stronger constraints can be obtained if NGC 4151 enters the flux 
level of the dim state for a longer period, as at A 2-10 kev - 10“®, 
Tph is close to the extent of the y-ray emitting region for thermal 
protons. 
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3.3. NGC 6814 

The y-ray signal from the direction of NGC 6814 with TS = 25.6 
was found bv lAckermann et al.l (l2012al) . However, the source is 
not reported in 3FGL. It may not have passed the TS > 25 crite¬ 
rion, assumed in 3FGL, as the above value is close to the limit. 
Our analysis, using data for 6.4 years, confirms the signal with 
TS ^ 32.3, and our best-fit parameters for a power-law spec- 
trum (see Table |2]) are a pproximately consistent with those of 
lAckermann et alJl2012al) . 

In Fig. [T^ we show the parameters of the intrinsic X-ray 
emission from spe ctral fits to Suzaku observations in 2011 
(IWalton et alJl20I^ . Interestingly, the source is located in the 
region for an inefficient (thermal synchrotron) source of seed 
photons. During this observation the source was in a relatively 
dim state, with the average ^wZ/f/BAT flux lower by a factor of 
~ 1.5 than the average during the Ferm//LAT observations, see 
Fig. Si. 

As noted above, we considered two values of M for NGC 
6814. For both, the measured Fi_ioGeV gives Ti-iocev exceed¬ 
ing the maximum value predicted for the hot flow (in models N 
with 6 — 10“^), which rules out an origin of the observed y- 
ray signal in an inner hot flow. For the higher value of M, NGC 
6814 has a similar /li_io cev to the y-ray-loud Seyfert 2 galaxies, 
a possible explanation is discussed in Sect.|5] 

3.4. NGC 5548 and 1C 4329a 

The likelihood analysis of the FermifLAI data does not reveal 
any significant signal around these two X-ray-bright Seyfert 
galaxies (in both TS < 1) and therefore we computed the ULs 
given in Table |2] 

For NGC 5548, we show in Fig. [T^ the parameters of the 
power-law fits from Brenneman et al. (2012) for seven Suzaku 
observations in 2007. The average SwiftlBA^ flux during the 6.4 
years of FennifLA^ observations is equal to the average flux dur¬ 
ing the Suzaku observations, and we used the intrinsic A 2 - 1 Q kev 
from the fit of Brenneman et al. (2012) for the average Suzaku 
spectrum. 

For IC 4329a, we show in Fig. [T^ the parameters of the 
power-law fits from Brenneman et al. (2014) for lower and 
higher flux states during Suzaku and NuSTAR observations in 
2012. We used /l 2 -i 0 kev larger by a factor 1.3 to account for the 
difference between the average SwiftfBAI flux for 6.4 years and 
the flux from measurement simultaneous with SuzakulNuSTAR 
observations, see Fig.|4j;. 

Both sources are observed above /lADAF.max but within the 
range of approximate applicability of our model. During one ob¬ 
servation, NGC 5548 is seen in the area of ADAF solutions with 
inefficient cooling. In the remaining observations both sources 
show interesting hints of modest hardening below or softening 
above /l 2 -iokeV - 0.002 with increasing /l 2 -i 0 kev- 

For both objects the ULs on the y-ray photon flux is lower 
by a factor of 2-3 times than for NGC 4151. In addition, both 
sources have a similar Ly/Lx ratio to NGC 4151 (for IC 4329a 
it is even lower). However, we see in Fig.[T]3 and Fig.[3]that they 
insignificantly constrain the hot flow models. 

3.5. NGC 7213 

Despite the large uncertainty in M, NGC 7213 gives interest¬ 
ing constraints on the hot flow model; crucially, it is observed 
at /l 2 -i 0 kev - 10“^, between the X-ray bright Seyfert 1 galax¬ 
ies discussed above and NGC 4258. The likelihood analysis of 


the FermifLAI data does not reveal any significant signal around 
NGC 7213 (TS < 1), and we computed the ULs given in Table 
|2] In Fig . [IK we show the pa rameters from a Suzaku observation 
in 2006 (iLobban et al.ll2010ll during which the average flux level 
was similar to the average during 6.4 years, see Fig.|4^. 

Figure[TJ) shows that the derived UL for Fi_io oev constrains 
the nonthermal energy content to rjp $ 0.1 for 6 - 10^^. Models 
with 6 - 0.5 predicts Ti-iocev below the Fermi limit, except 
for a - Q and rjp - \. However, the X-ray spectra predicted for 
6 - 0.5 appear to be qualitatively inconsistent with the spectrum 
measured in NGC 7213. The spectra at /l 2 -iokeV - 10^^ (the 
dotted line for a - 0.95 and the overlying dashed line for a = 0 
in Fig. 12 };) deviate significantly from a power-law shape, whereas 
the Suzaku spectru m between 0.6 and 5 0 keV is described by a 
simple power-law (iLobban et al.ll201^ . The UL for F' 0 . 2-1 Gev 
rules out a high (5 flow around a high-a black hole. 

3.6. NGC 4258 

The precisely determined parameters of this AGN made it an es- 
sential object for t he development of the hot flow models (e.g., 
iLasota et al.lll996h . However, it has not been included in pre¬ 
vious analyses of the Fermi/LAI data. At the position of NGC 
4258 we found a weak y-ray excess above the background, with 
TS = 9.4. A likelihood analysis of the signal with gtlike gives 
F = 2.5 + 0.2, F ^ (2.8 + 1.2) x 10-“^ ph s^' cm^^. Despite 
the weak residual, the derived ULs on the photon flux yield the 
tightest constraint on /li_io GeV and /I 0 . 2-1 GeV of all considered 
objects, see Table |2l 

In Fig. dH we sho w parameters of the analysis reported by 
iRevnolds et al.l (l2009h of the Suzaku observation in 2006, who 
also noted that during this observation NGC 4258 increased its 
intrinsic X-ray luminosity by a factor of ~ 2 relative to the aver¬ 
age from Swift/BAI measurements. We used /l 2 -iokeV from the 
Suzaku fit reduced by a factor of 2. 

At /l 2 -iokeV - 10'^, the 0.2-1 GeV emission is weakly af¬ 
fected by yy absorption, and we can observe radiation from the 
ergospheric region (see Sect. 2). However, as Fig. [3] shows, the 
strongest expected signal is only at the level of the Fermi UL. 

For 6 — 10^^, model N predicts a flux strongly exceeding the 
derived UL for Fi_io gcVj ruling out rjp ~ For 5 - 0.5, the pre¬ 
dicted y-ray flux is an order of magnitude below the UL value. 
However, the corresponding X-ray spectrum (bottom dashed line 
in Fig. 12 };) seems to rule out this (large-b) case if the X-ray radi¬ 
ation observed from NGC 4258 is produced by thermal Comp- 
tonization. 

3.7. Circinus, NGC 1068, and NGC 4945 

For the Seyfert 2 galaxies we used the parameters re¬ 
ported in 3FGL, which are approximately consistent (F i mGsv 
within 30%) with previous analyses (iLenain et aTI 1201 Ot 
lAckermann et al.l[2012bt iHavashida et al.ll2013h. 

Of these three AGNs, only NGC 4945 shows rapid variabil¬ 
ity in the hard X-rays, which indicates that it is a transmission- 
dominated source. We can therefore d irectly probe the nuc lear 
emission. In Fig.[T^ we show fits from iPuccetti et all (1201 4ll for 
NuSTAR observations in 2013. The average SwiftfBAI flux dur¬ 
ing the four years of Fermi observations used in 3FGL agrees 
with the average flux during the NuSTAR observations. 

In NGC 1068, the direct X-ray emission of is completely 
obscured along our line of sight and we only see the reflected 
component. Then, the assessed intrinsic emission is strongly 
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m odel dependent. We used A 2-10 kev corresponding to the best-fit 
of iBauer et al.l (l2014t) . This assessment, based on the observed 
reflected component, gives the level of emission from the ac¬ 
tive nucleus averaged over a long time, possibly over hundreds 
of years, because a significant fracti on of reflection arises at a 
~ 100 pc scale (see lBauer et al.l2014l) . Similarly, the X-ray emis- 
sion from Circinus is reflection dominated; we used results from 
lArevalo et al.l (l2014h . NGC 1068 and Circinus are not shown in 
Fig. [1^ because the estimates for their intrinsic Fx are model 
dependent. 

Figure [TJ) shows that all three y-ray-loud Seyfert 2 galaxies 
agree remarkably well after scaling by their central black hole 
masses, with A 2 -wkeV - 0.02 and Ti-iocev - 3 x 10“^ in all 
three. We also note that their Ti -10 Gev/T 2 -io keV ratio does not 
exceed the limit on the luminosity ratio in other AGNs, except 
for that corresponding to the UL in IC 4329a for F = 2.7. 

3.8. Centaurus A 

The interpretation of the high-energy spectrum of Cen A has 
some open questions. Its X-ray emission may contain contribu¬ 
tion from both the jet and the accretion flow. T he latter origin is 
suppor ted by the thermal-like cutoff, claimed in iBeckmann et al.] 
(1201 ih . and by the location in the area of the /I 2-10 keV-Fx plane 
occupied by Seyfert galaxies. In Fig. [TK we show the parameter s 
from three Suzaku obs ervations in 2009 (iFukazawa et al.ll201 ih . 
iFukazawa et al.l (l20I ih also noted hints for a harder power-law 
component, with Fx < 1.6, weakly contributing below 100 keV, 
which may represent emission from the jet. The main contribu¬ 
tion to the r-r ay emission of Cen A most likely comes from the 
jet. However, ISahakvan et al.l (1201 3h found evidence for a sec¬ 
ond component, an order of magnitude less luminous than the 
main y-ray component. This weaker component has a hard spec¬ 
trum and dominates above ~ 4 GeV. Assuming that this com¬ 
ponent extends down to 1 GeV with the same slope, we find 
di-ioGeV - 1.7 X 10~®. Based on a comparison with predictions 
of model Ho.i , this component can be emitted from a flow with 
6 - 10“^ and rjp =^Q.l - 0.2 (slightly dependent on plasma mag¬ 
netization). 

The SwiftfBKT flux level during the 2009 Suzaku observa¬ 
tions is equal to the average flux d uring the FermifLKf d ata- 
taking period used for the analysis in ISahakvan et al.l (l2013h . 

4. Nonthermal electrons 

Here we consider effects related with nonthermal electrons that 
are produced by various mechanisms, specifically, their contri¬ 
bution to the observed y-ray luminosity and to the optical depth. 

4.1. P/on decay 

If nonthermal protons are present, electrons with y ^ 10^ may be 
produced by charged pion decay. The nonthermal Compton radi¬ 
ation of these electrons may extend to the LAT energy range. The 
power injected in decay electrons is similar to the luminosity 
in 7r° decay y-rays. However, electrons with y ^ 10^ lose most of 
their energy in synchrotron radiation, therefore their nonthermal 
Compton emission in the y-ray range is always much weaker 
than the decay component. 

The density of -decay determined by the pair equilib¬ 
rium condition (cf. iNiedzwiecki et al.ll2015ll . is at least two or¬ 
ders of magnitude lower than the density of the ionization elec¬ 
trons. The production rate of is similar to that of there¬ 


fore pairs produced by absorption of Tr^’-decay y-rays give at 
most a similar contribution to the total density as those from n"^- 
decay. Production of further generations of pairs is unlikely 
for the luminosities described by our model, see below. Overall, 
the products of hadronic processes contribute negligibly to the 
optical depth. 

If a strong yy absorption affects the ;7r°-decay y-rays, the syn¬ 
chrotron radiation of thus produced e*- pairs increases the seed 
photon input, again at most by a factor of ~ 2. We neglected this 
effect in computations presented in this work; we checked that 
this twice higher seed photon input yields an electron tempera¬ 
ture that is reduced by up to 10% and only insignificantly affects 
our results. 

4.2. Direct acceleration 

We briefly discuss constraints on the acceleration of electrons 
resulting from the FermifLAI ULs. For the results presented in 
this section we do not use the simplified description (i.e., with a 
monoenergetic injection) of the accelerated electrons, and we as¬ 
sume that they are accelerated to the power-law distribution with 
the acceleration index s^. As above, the amount of nonthermal 
electrons is parametrized by rj^. Specifically, the power used for 
the nonthermal acceleration is Qacc = ?7eb2diss, where Qdiss is the 
total power dissipated in the flow. 

Conditions in the inner parts of hot flows surrounding super- 
massive black holes allow for accelerati on of electrons t o GeV 
energies; for instance, using Eq. (7) from IZdziarski et ^ (l2009ll . 
we found the maximum Lorentz factors to be limited by the syn¬ 
chrotron losses of ymax ~ 10^-10® for the typical magnetic field 
strength of B ~ (10^ -10"*) G. The synchrotron emission of these 
electrons extends up to ~ 10 MeV, so they can produce photons 
in the LAT energy range only by Compton scattering. We im¬ 
plemented the description of nonthermal Compton emission in 
our MC code and performed simulations with ymax = 2 x 10^, 
for which the nonthermal Compton spectrum extends up to ~ 10 
GeV. 

The y-ray luminosity produced by the nonthermal Compton 
scattering primarily depends on and s^. We found that for 
ie ^ 2.4 and rje = 1, the values of /I 0 . 2-1 GeV and /li_ioGeV for 
the nonthermal Compton are below any UL derived from the 
LAT data for objects with /l 2 -iokeV ^ 10“^. Lor ^ 2.4, the 
y-ray Compton component for - 1 exceeds the ULs for NGC 
4258, NGC 7213 and NGC 4151, then, for such ie the fraction of 
accretion power used for nonthermal acceleration is constrained. 
We found for instance that for ie - 2 it is limited to ^ 0.052diss 
(corresponding to b = 0.5 and S 0.1). 

In our computations we assumed that synchrotron radiation 
emitted below the self-absorption frequency is thermalized, as 
described in several works on the synchrotron boiler mecha¬ 
nism (e.g.- lMalzac & Belmonlll2()()^ . meaning that the effective 
heating of thermal electrons is higher than (1 - ?7e)b2diss- The 
power in the self-absorbed radiation is ^ 0.22acc for Se ~ 2. 
Lor Se = 2.6, it is ~ 0.82acc- Then, if 6 is large, the heating of 
thermal electrons is always strong, with the synchrotron boiler 
dominating for ~ 1. This also implies that our conclusions re¬ 
garding the allowed values of are weakly dependent on other 
parameters of the flow. Although Qdiss shows a significant de¬ 
pendence on a, for example, both the nonthermal Compton Ly 
and the thermal Compton Lx scale roughly linearly with Qdiss 
(the latter either through the synchrotron boiler or direct MHD 
heating), making the LyjLx ratio weakly dependent on the other 
parameters. 
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4.3. yy absorption 

Secondary pairs, produced by yy absorption, in general affect 
radiative properties of flows in the considered range of A 2 - 1 Q kev 
only insignificantly. The first generation of pairs contributes at 
most a few per cent to the Thomson optical depth. Copious cre¬ 
ation of secondary pairs through electromagnetic cascades re¬ 
quires (1) the Compton process to dominate over synchrotron, 
and (2) Compton photons produced above the threshold for pair 
production. Typically, at least one of these conditions is not 
fulfilled. In the innermost parts of the flow, Compton energy 
losses are comparable to synchrotron losses only for electrons 
with Lorentz factors not exceeding a few hundred. At higher en¬ 
ergies, synchrotron losses strongly dominate as a result of the 
Klein-Nishina decline in the Compton scattering rate, hence, 
the development of a cascade is suppressed. In turn, the lower 
energy electrons, for which the more efficient Compton losses 
are mostly due to interactions with X-ray photons from thermal 
Comptonization, mainly produce photons with £ ^ 10 MeV. 
However, at /i 2 -iokeV ^ 10“^, the flows are mostly transparent 
to such photons. Therefore, we do not expect a significant con¬ 
tribution of secondary pairs to the Thomson depth of the flow. 
We note, however, that at /i 2 -iokeV ^ 10“^ the innermost region 
should be opaque to photons with E k 1 MeV, and in this regime 
the secondary pairs may then contribute considerably to the op¬ 
tical depth. 

We also note that spatially exte nded pair cascades, such as 
discussed by IZdziarski et alJ (l2009h . for instance, are unlikely 
to strengthen the observed y-ray flux. Such extended cascades 
can only develop in the region with Tyj ~ 1. In a hot flow, Tjy 
increases quickly toward the center, see Fig. 7 in N13, there¬ 
fore the free path of y-rays produced within rph is typically very 
short. Furthermore, as noted above, the absorbed energy is lost to 
synchrotron radiation rather than used for cascade development. 


5. Discussion 

We first briefly comment on the X-ray spectral properties of 
the considered AGNs with /I 2-10 ^ 10^^. At these luminosities, 
the seed photons for Comptonization must be produced inter¬ 
nally to the flow, and nonthermal synchrotron radiation is an 
almost obvious process providing the amount of seed photons 
needed to explain the A GN data (iPoutanen & Veledinal 120141 : 
iNiedzwiecki et alJl2015h . If 6 is small, the nonthermal electrons 
are provided by the decay of charged pions, unless both a and 
P are small and protons are thermal, see Sect. 12.11 We empha¬ 
size that the weakness of y-ray signals related with decay 
does not contradict a copious injection of nonthermal electrons 
by decay because the most efficient pion production occurs in 
the region opaque to y-rays. If S is large, the decay does not 
provide a sufficient amount of nonthermal electrons because the 
flow density is too low (here inefficient production is strictly 
related with very low level of y-ray activity). Then, direct ac¬ 
celeration must take place if 6 is large. These models (i.e., with 
either a significant production or T/e ^ 0.1) All the upper gray 
region in Fig.[T^. Typical Seyfert spectra clearly agree with these 
solutions. 

Solutions with seed photons from thermal synchrotron alone 
(including models with 6 - 0.5 and - 0, as well as model T 
with 5 - 10“^, a - 0 and >6=1) fill the lower gray region. Ob¬ 
viously, intermediate cases are also possible, for instance, with 
5 - 0.5 and //e ^ 0.01. Seyferts are occasionally observed with 
parameters corresponding to inefficient cooling (NGC 6814 and 
one of NGC 5548 data points). These harder-than-typical spectra 


may result from a decline of nonthermal electron injection. For 
large 6 it may correspond to the decrease of rje, and for small 5 
to the decrease of >6 (again, for small a and thermal protons). 

Observational assessment of nonthermal acceleration pro¬ 
cesses in astrophysical plasmas is an essential question. Hot 
flows in which ultra-relativistic acceleration of protons takes 
place can be considered a sour ce of neutrinos detected by Ice- 
Cube (e.g.. lKimura et al]|2015h . Accretion flows should also be 
strong acceleration sites according to hybrid thermal and non¬ 
thermal Comptonization models (e.g., Vurm & Poutanenll2009l : 
iMalzac & Belmoti^ 120091 : IVeledina et ^ 201 Ih that are widely 
applied to explain X-ray spectra of accreting black holes; here 
at least ~ 50% of the energy must be provided to electrons by 
nonthermal processes if thermalization of electrons via Coulomb 
collision and synchrot ron self-absorption is taken into account 
(iDel Santo et alj|2013h . Interestingly, these models are typically 
fitted with ^ 2.4. As discussed in Sect. 14.21 the y-ray lumi¬ 
nosities from nonthermal Compton emission of flows with such 
Se are below the LAT ULs deriyed in Sect.[3]for Seyfert galaxies, 
eyen if — \. 

Obseryational eyidence of the presence of non-thermal par- 
ticles in accretion flows includes detection of MeV tails (e.g., 
iDroulans et al.ll2()T(il : IZdzia rski et al' 2012h a nd patterns of op¬ 
tical/infrared eyolution (IPoutanen et al.ll2014t) . both obseryed in 
stellar black-hole systems. Howeyer, these features are produced 
by nonthermal electrons, which may come from pion decay, and 
in principle, they do not require direct acceleration. 

Obseryations in the GeV range proyide the means to directly 
measure the efficiency of proton acceleration in hot flows. The 
obtained ULs for £ 1-10 Cev in NGC 4258, NGC 7213 and NGC 
4151 constrain the fraction of accretion power used for relatiyis- 
tic acceleration of protons to at most 77p ^ 0.1 (if 5 ^ 0.1). Re¬ 
markably, these three AGNs allow probing yarious parts of the 
flow because at /i 2 -iokey ~ 10“^, 10^“* and 10“^, radiation pro¬ 
duced at a few, ten, and seyeral tens of £g, respectiyely, domi¬ 
nates the obseryed 1-10 GeV flux. The component dominating 
the core emission of Centaurus A aboye ~ 4 GeV may come 
from a hot flow in which ^ 0.1 - 0.2. In this scenario, its 
T/p should be higher at least by a factor of ~ 2 than in NGC 
7213 (see Fig.[T}5), which has similar 42-10 keV- This might reflect 
differences between flows powering radio-loud and radio-quiet 
AGNs. We also note that in this hot flow scenario for Centaurus 
A, the obseryed y-rays (aboye 4 GeV) proyide information about 
the region close to the jet formation sit^B 

We now briefly discuss hard spectral states of black hole 
binaries in which hot accretion flows, similar to those power¬ 
ing AGNs, are most likely present. In the innermost parts of 
flows powering both classes of objects the same physical pro¬ 
cesses are expected. This yiew is supported by the similarity be¬ 
tween the X-ray spectra of S eyferts and black hol e binaries in 
the hard spectral states (e.g.. IZdziarski et^ll996ll . with slight 
differences simply explained b y scaling of seed photon energies 
with the black hole mass (cf. INiedzwiecki et al.ll2015h . These 
differences do not affect the y-ray emission, and the relation be¬ 
tween /I 2-10 keV and 4o.2-i CeV or 4i_io cey shown in Figs[TH3]is 
probably also releyant for flows around stellar-mass black holes. 

An estimate of the y-ray actiyity leyel, which can be com¬ 
pared with the hot flow mo del, has only be e n obta ined for the 
hard state of Cyg X-1 by iMalyshey et al.l (l2013h . Figure [TJr 
shows that it giyes a slightly tighter constraint than Seyfert galax¬ 
ies, limiting r/p to at most a few per cent for 6 - 10“^, and at 


^ according to the radio imaging of M 87, e.g. lDoeleman et ^ ( l2012h . 
and assuming that the sizes of such sites are common in AGNs 
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most ~ 10% for 6 - 0.5 and a - 0. A significant nonthermal 
proton component, with rjp ~ 0.5, is only allowed if 5 ~ 0.5 
and a is high. This finding places interesting constraints on hy¬ 
brid models (see above) if protons are accelerated with at least a 
similar efficiency to electrons (as expected in many acceleration 
scenarios). It may be also interesting to note that similar values 
of /li_io Gev correspond to the Cyg X-1 hard state (^ 6 x 10“^) 
and our gtlike fits to signals at the positions of NGC 4151 
(^ 5 X 10-^) and NGC 4258 (^ 3 x 10^’), although the signifi¬ 
cance of the y-ray signal is rather weak (TS = 15.6) in Cyg X-1 
and marginal in NGC 4258 and NGC 4151. 

Finally, we comment on the y-ray emission of NGC 4945, 
NGC 1068, and Circinus. All three exhibit both an AGN and 
starburst activities. NGC 4945 and NGC 1068 approximately 
follow some multiwa velength correlations esta blished in other 
star-forming galaxies (lAckermann et al.ll2012bl) . which suggests 
that their y-ray emission may be related to the latter activity. 
However, Circinus does not agree with these correlations, fur¬ 
thermore, its Ly exceeds the calorimetric limit (where all cosmic 
ray s produced by supernov ae interact wit h interstellar matter) , 
see iHavashida et al.l (l2013h . Furthermore, iLenain et al.l (1201 (ih 
noted that Ly measured in NGC 1068 exceeds by a factor of 
~ 10 the level expected for the total gas mass and supernova rate 
estimated in this galaxy, which led the m to suggest that the y- 
rays may be produced in a jet. Similarlv. lYoast-Hull et al.l(l2014h 
found that their detailed starburst model underestimates the ob¬ 
served y-ray flux and overestimates the radio flux for NGC 1068, 
and they argued that these problems would be resolved if the 
AGN were the primary source of y-rays. This shows that the ori¬ 
gin of y-rays still seems rather unclear. 

Seyfert galaxies show a strong correlation between 
nuclear star-formation and t he AGN luminosity (e.g., 
iDiamond-Stanic & Riekel l2012h . Masses of supermassive 
black holes in Circinus, NGC 1068, and NGC 4945 are pre¬ 
cisely determined by H 2 O megamaser measurements. Their 
values imply an intriguing similarity of the values of both 
T 2 -iokeV and Ti-iocev in all three objects. The question then 
arises whether their y-ray emission might be related with 
accretion processes in their active nuclei. For the intrinsic 
42-iokeV - 0.02, the geometry of accretion flow is very uncer¬ 
tain. It cannot be studied through Fe Ka line distortions because 
the X-ray spectra of all three AGNs are completely absorbed 
below 10 keV. Some analogies with black hole binaries, in 
which (1) lu minous hard states ha ve been observed at such 
42-10keV (e.g.. IZdziarski et ani2004h . (2) softening of intrinsic 
X-ray slopes above /l 2 -iokeV - 10“^ is observed that closely 
resembles t he beh avior seen in NGC 4945 (see, e.g.. Fig. 1 in 
lYang et ^ (1201 5l) I. and (3) strong thermal components are seen 
at /I 2-10 keV ~ 0.01, indicating that t he cold disk extends down to 
rtr ^ 10 te.g.. lCabanac et alll2009l) suggest that inner hot flows 
may be present in these AGNs, but only within the innermost 
several Rg. This value of ~ 10 is tentatively shown in Fig. 
[Di. For such a small, strongly energized (the main release of 
gravitational energy occurs in this region) and relatively dense 
flow, we cannot exclude the possibility that a sufficient y-ray 
flux is produced in an outer layer and escapes avoiding yy 
absorption (taking into account a strong local anisotropy of the 
X-ray radiation field, see Sect. 2). Any stronger conclusions 
for this scenario would be premature; a self-consistent model, 
compatible with the X-ray data, needs to be developed before a 
quantitative assessment of the y-ray flux can be made. 

We do not expect direct y-ray signals from central parts of 
accretion flows in black hole binaries even if they have the same 
geometry at this A 2-10 keV- Any emission in the GeV range would 


be absorbed in interaction wi th soft X-ray em ission of a cold 
disk extending to small r, see iBednarel^ (Il993h for quantitative 
estimates. In AGNs, a UV-emitting cold disk only attenuates y- 
rays at ~ 100 GeV. 

Another possibility involves escape of protons accelerated in 
the flow and their interaction with circumnuclear matter. About 
0.1% of the accretion power would need to be carried away by 
escaping protons to explain the y-ray measurements. This would 
also allow explaining the y-ray signal in NGC 6814 by the same 
emission mechanism, assuming that (1) its supermassive black 
hole has M ^ 2 x 10^Mq (the higher of the two values consid¬ 
ered above), for which its /li_ioGeV is similar to those of NGC 
4549, Circinus, and NGC 1068; and (2) protons were injected 
into the circumnuclear region during an earlier phase of higher 
(by almost 2 orders of magnitude relative to current) activity 
of NGC 6814, in th e scenario similar to that proposed for the 
Galactic center (e.g.. lChernvakova et al.ll20Tlh . In this scenario, 
an energy-dependent diffusive loss of protons may explain the 
slightly softer y-ray spectrum of NGC 6814. 

6. Summary 

Hot flows are supported by the proton pressure, and the produc¬ 
tion of pions in collisions of the energetic protons should be their 
generic property. If this is true, a lack of detection of a related 
y-ray signal could be regarded as an argument against this class 
of models. However, the expected flux depends on several pa¬ 
rameters, and in the extreme case of a strongly magnetized flow 
around a nonrotating black hole, the conditions for pion produc¬ 
tion are not achieved if protons are thermal. 

We compared the prediction of the model with several well- 
studied AGNs, for which the available data allowed a robust de¬ 
termination of the nuclear luminosity scaled by the Eddington 
value. At the positions of NGC 4151 and NGC 4258, we found 
weak residuals, whose nature (background fluctuation or an ac¬ 
tual emission) could not be assessed because of their low statis¬ 
tical significance. 

The highest y-ray luminosities are predicted if electrons are 
heated mostly by Coulomb interactions. If additionally most of 
the accretion power is used for the relativistic acceleration of a 
small fraction of protons, the highest level of /li_ioGeV - 10“^ 
is expected for /l 2 -iokeV between ~ 10“^ and 10“^. This exceeds 
the Fermi limit on /li_io cev for NGC 4258, NGC 7213, and NGC 
4151 at least by a factor of several. Thus, the Fermi ULs provide 
an interesting constraint on the MHD processes that convert the 
accretion power into the kinetic energy of protons: A rather uni¬ 
form heating of all protons would need to take place. Given this 
constraint, we used Ao. 2-1 GeV predicted for thermal distribution 
of protons to rule out a flow with a high (equipartition) value of 
yS around a high-a black hole for NGC 7213 and and NGC 4151; 
for the latter this result is model dependent. 

If a large part of the accretion power is used for a direct heat¬ 
ing or acceleration of electrons, the expected y-ray fluxes are 
below the Fermi ULs, except for NGC 4151 and NGC 7213, for 
which T]p ^ 0.5 was ruled out with a small a (yielding some con¬ 
straints on the acceleration efficiency, which may be relevant for 
hybrid models). However, this (large 6) class of models appears 
to be disfavored by the predicted properties of X-ray spectra, 
which significantly deviate from a power law at A 2-10 keV ^ 10^^ 
and are described by (too) small t at larger A 2 -W keV- 

The y-ray-loud Seyfert 2 galaxies radiate at a much higher 
Eddington ratio than other nearby AGNs, possibly as a result of 
enhanced fueling that is related with starburst activity in their 
nuclear regions. At this /l 2 -iokeVj the transition between the hot 
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flow and the cold disk possibly takes place within the central 
107?g. Whether it creates conditions suitable for the escape of 
y-rays or relativistic protons remains an open question. 
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